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ABSTRACT 

Recent data indicates that the cosmic UV emissivity decreased with decreasing redshift 
z near the end of reionization. Lacking evidence for very massive early stars, this could 
signal a decline with time in the mass-averaged escape fraction of ionizing radiation 
from galaxies (f esc ) at z 6. We calculate the evolution of ionization fronts in dark 
matter halos which host gas in hydrostatic equilibrium at its cooling temperature 
floor (T w 10 4 K for atomic hydrogen). We find a high escape fraction only for the 
lowest mass halos (with M < 1O 8 - 7 M at (1 + z) = 10) provided their star formation 
efficiency /* <L 10~ 3 . Since the low-mass galaxy population is depleted by radiative 
feedback, we find that indeed (/ esc ) decreases with time during reionization. 
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1 INTRODUCTION 

The first generation of galaxies is expected to have reion- 
ized cosmic hydrogen by a redshift z « 7 ( |Loeb fc Furlan-] 
|etto|2013 1. One of the most important unknown parameters 
regulating the way in which reionization proceeds is the frac- 
tion of ionizing photons, / esc , that escape outside the virial 
radius of the galaxies in which they were produced and into 
the surrounding intergalactic medium (IGM). The growth 
of ionized regions with cosmic time depends on the average 
value of this parameter over viewing angle per galaxy and 



over the galaxy population (Wyithe & Loeb 2003 Trac & 
|Cen|2007| |Bouwens et al.|2012| |Finkelstein et al.||2012| ). 

Observations of galaxies at 2 = 1 — 3 indicate a broad set 
of escape fraction values ranging from a few percent to tens 



of percent ( 


Steidel et al.||2001 


Fernandez-Soto et al.||2003 


Shapley et al.||2006| |Inoue et al.||2006. 5 


jiana et al.||2007 
outsia et al.|2011 


Giallongo et al.|2008||Iwata et al.|2009||B 


Grazian et al. 2011 


Vanzella et al. 2012a b 


I. This potentially 



reflects strong variations of the escape fraction with viewing 
angle and evolutionary time. 

A calculation of the escape fraction of ionizing radiation 
from first principles is difficult, as it depends on the spatial 
distribution of ionizing sources relative to the neutral hydro- 
gen in galaxies and is therefore sensitive to the small scale 
clustering of young stars and the interstellar medium. Due to 
their higher density, high-redshift galactic disks are expected 



to have small escape fractions during reionization 1 


Wood & 


Loeb|2000 ; Gncdin et al. 2008 , Razoumov & Sommer-Larsen 


2006 Wise & Cen|2009 


1. Even for local disk galaxies, 


Dove 



|et al.| ( |2000[ ) concluded that f esc should not exceed a few 
percent in most cases, as a result of efficient radiation trap- 
ping by the shells of the expanding superbubbles around OB 
star associations. 

Since the mean density of the Universe (and hence 
galaxies) scales as (1 + z) 3 , theoretical calculations tend 
to conclude that f esc should decrease with increasing red- 
shift. This expectation appears to be in conflict with recent 



reionization models and data ( Bolton & Haehnelt|2007| |Sal- 


vaterra et al. 


2011||Mitra et al.J2012a|b||Kuhlen & Faucher- 


Giguere|2012 


Finkelstein et al.|2012 1 that point towards the 



need for an increasing UV emissivity towards high redshifts. 
This can be achieved in two ways. The first is to postu- 
late an increasingly top-heavy initial mass function (IMF) 
of stars that would increase the number of photons emitted 
per baryon incorporated in stars. Although possible, this ex- 
planation seems to be disfavored by the lack of observational 
evidence for an early population of substantially more mas- 



sive stars (Cayrel et al. 2004 Caffau et al. 20111 and by 



recent numerical simulations ( Greif et al.|2012 1. A more ap- 



pealing solution is represented by a possible increase in the 
average escape fraction (/ eS c} during reionization, recently 
suggested by Mitra et al. (2012b I. Here we explore a novel 



physical explanation for this unexpected trend. 

We show that the lowest-mass galaxies near or below 
(so-called mini-halos) the hydrogen cooling threshold have 
high escape fractions but their contribution to the UV pro- 
duction rate decreases with cosmic time due to large-scale 
radiative feedback processes during reionization that either 
photo-heat or sterilize them by dissociating their H2. How- 
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ever, internal feedback is of foremost importance as well: as 
these small systems achieved large f esc by rapidly (within 
~ 10 4 yr) ionizing their entire interstellar medium, this gas 
will become loosely bound to the host galaxy. Under these 
conditions it is difficult to sustain a continuous mode of star 
formation (SF), particularly because following the death of 
massive stars powerful supernova explosions will clear the 
gas out of the potential well before a stable disk structure 
is established. More massive galaxies will form a disk with 
the net result of trapping most of their ionizing photons. 

If high-redshift dwarf galaxies form stars predominantly 
in episodic bursts, as suggested by the high incidence of 
galaxy mergers at increasing redshift ( Munoz fc Loeb|2011[ ) , 
the conventional argument that their star formation effi- 
ciency is suppressed by supernova (SN) feedback may not 
be valid. In particular, SN feedback is inhibited if the dura- 
tion of the starburst is shorter than ~ 10 Myr, the lifetime 
of SN progenitors. In this regime, the global radiative feed- 
back sets the amount of gas initially available for making 
stars and hence the overall efficiency of star formation, 
In our calculations, we will calibrate the value of /* based 
on recent measurements of the cosmic stellar mass density 
at high redshifts ( Gonzalez et al.|2011 1. 

For an instanteneous starburst, / esc can be large only if 
an escape route for ionizing photons is opened within a few 
Myr, prior to the death of the massive stars that produce 
these photons. This process cannot be mediated by super- 
nova explosions, which occur after the emission of ionizing 
photons has already started to decline. The key question is 
therefore whether the surrounding blanket of absorbing hy- 
drogen atoms can be ionized before the massive stars in the 
starburst end their life. We address this question within the 
context of a simple model in which we populate dark mat- 
ter halos with primordial gas that cooled to the temperature 
floor of atomic hydrogen, ~ 10 4 K. 

Our basic point is simple: massive galaxies exist at all 
redshifts and the escape fraction of ionizing radiation from 
them is calculated and observed to be small; however, the 
lowest-mass galaxies (T V i r ,$ 10 4 K) exist only before and 
during reionization and have a high / esc . Our assumption of 
a short-lived star formation phase without SN feedback is 
self-consistent, as we will show that photoionization heating 
of the gas suppresses star formation in dwarf galaxies after 
a period of time much shorter than the lifetime of massive 
stars. 

In §2 we describe the details of our calculation and in 
§3 we present our numerical results. Finally, we discuss our 
main conclusions in §4. Throughout the paper, we adopt the 
WMAP7 set of cosmological parameters with Q. m = 0.27, 
Q.A = 1-H m , n b h 2 = 0.023, h = 0.71, a 8 = 0.81, n s = 0.97, 



and dn s /dlnfe = (La rson et al.|20lT l 



2 METHOD OF CALCULATION 

We derive the escape faction of the ionizing photons, f esc , 
based on a set of simplifying assumptions. We consider a 
spherical dark matter halo of mass, M , virializing at redshift 
z and characterized by an internal density (spherically av- 



Table 1. Values of central number density, no (in units of cm -3 ) 
for dark matter halos of different total mass, M, at z = 9 and for 
different values of the gas temperature, T, expressed in terms of 
the virial temperature T V % T through the parameter V 2 = T v i r /T. 



log 10 T„ ir [K] 


4.31 


4.98 


5.65 


6.32 


M [M ] 


10 8 


10 9 


10 10 


10 11 


V 2 = 1.0 


20.6 


23.9 


36.0 


167.8 


V 2 = 1.2 


46.2 


54.5 


87.3 


282.9 


T = 2 x 10 4 K 


24.3 


1.9 x 10 4 


4.3 x 10 6 





|et~aLl|T997| ), NFW) form, 
p(r) = 



pcSc 



(1) 



cx(l + cx) 2 ' 

where x = r/r V i r , r V i r is the virial radius of the system, c is 
the halo concentration parameter, 8 C = (18n 2 / 3) c 3 / F(c) is 
a characteristic overdensity, and 

F(i)=ln(l+t)- r ^. (2) 



For the concentration parameter we use the results of Prada 
et al. (20121, extrapolating them towards lower masses 



and higher redshift when necessary. The circular velocity, 
v 2 (r) — GM(r)/r, can be expressed as 

«c = K 2 ^, (3) 

xb (c) 

with V 2 = GM/r V i r . Baryons fall into the dark matter halo 
potential well are shock-heated to the virial temperature 
T vir = {pm p /2k)V 2 , and then relax to a hydrostatic config- 
uration ( Makino et al.|1998 l, 

V 2 2 2 
In p(r) = In p - — [v e (r = 0) - v e (r) . 

Here, v e is the halo escape velocity from the halo, 
r 2 F(cx) + cx/(l + cx) 



(4) 



v 2 (r) 



2v; 



xF(c) 



(5) 



p — 1.22 is the mean molecular weight of a neutral primor- 
dial H-He gas; ks is the Boltzmann constant and m p the 
proton mass. We have also inserted in equation Q the extra 
factor V 2 = V 2 /c 2 s = T vir /T > 1, where c 2 s = 2k B T/pm v is 
the effective sound speed at a generic temperature T / T vir . 
The standard case in which T = T V i r is obtained by setting 
V = 1. By manipulating equation Q we then get the final 
expression for the gas density profile as a function of radius, 



p(x) = poe 7 (1 + cxy /cx = pm v nor)(x)\ (6) 

where we have defined 7 2 = 2c/F(c), and the ancillary func- 
tion t](x) = e~ 7 v (1 + cx) 1 v ' cx . The central density, po, 
is obtained by requiring that the total gas mass in the halo 
is equal to the cosmological value, i.e. Mb = (S7f,/f2 m )M = 
fbM. This procedure yields 

po = (18n 2 /3)f b c 3 e^ V * 



dt(l + t) 



j' ! V' 1 /t^2 



(7) 



eraged) profile following the Navarro- Frenk- White ( Navarro 



Table 1 provides the values of the central number density, 
no, for different halo masses and V values at a fixed red- 
shift, z — 9. As expected, no increases in larger masses even 
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in the case V = 1 in which the gas is at the virial tempera- 
ture of the halo. This reflects the strong dependence of no on 
the concentration parameter, which increases from c w 4 for 
M = 10 s M to c w 6.5 for M = 10 11 M Q . This increasing 
trend is amplified if the gas is cooler than the virial temper- 
ature (V > 1) or if it is thermostated at a fixed temperature 
(such as T — 2 x 10 4 K in Table 1), leading to unrealis- 
tic densities for the most massive halos with M = 10 1 Mq 
(which represent rare 3a fluctuations of the density field at 
z = 9). For massive halos we obtain artificially high den- 
sities because we ignore angular momentum. Once the gas 
condenses to a sufficiently small scale, its rotation will halt 
contraction and it will settle to a disk. Disks were consid- 



ered by Wood fc Loeb (2000) and Dove et aL ( 2000 1 who 



independently concluded that the escape fraction was neg- 
ligible for a smooth gas distribution in high-redshift disks 
and also local spirals. We therefore argue that once the gas 
contracts enough to make a disk, its density will not increase 
by as much as our spherical model predicts but UV photons 
will not be able to escape from it anyway, based on earlier 
studies. 



2.1 Ionization front evolution 

Having specified the gas density distribution inside halos, 
we may now derive the evolution of an ionization front (IF) 
driven by the emission of photons with energy > 13.6eV 
from a burst of star formation at the halo center. We as- 
sume that a fraction /* of the available gas, fbM, of a given 
halo is instantaneously converted into stars distributed in 
mass according to a Salpeter IMF in the range (mi om , m up ) 
= (1M Q , 100 M Q ) and with absolute metallicity Z = 10~ 3 . 
The ionizing photon production rate, Q(t), by the stellar 
cluster can be computed exactly fr om population synthesi s 



models: we use here Starburst99 Leitherer et al. 



(1999) 



The time dependence of the production rate of ionizing pho- 
tons rate under these conditions is, 



Q(t) = 



Qo 

i + (t/t y 



(8) 



with (Qo,to) = (10 47 s _1 M Q -\ 10 6 - 6 yr). Equation Q illus- 
trates the important point that after ~ 40 Myr, the produc- 
tion rate of ionizing photons rapidly drops as a result of the 
death of short-lived massive stars. This has important im- 
plications for f esc as discussed hereafter. Note that equation 
|8| implies that the number of ionizing photons emitted per 
baryon incorporated into stars is N-, « 0.5Qtom p / Mq = 
5 x 10 3 . 

The time evolution of the IF radius, rr, is described 
by an ordinary differential equation that expresses a de- 
tailed balance between the ionization and recombination 
rates within the volume enclosed by the IF, 



dri 



47m jjr^ 



Q(t) 



Itt 



3 2 (2) 

rjn H a 



(9) 



where Tin = 0.92n is the hydrogen density for a primor- 
dial gas and q (2) = 2.6 x 10~ 13 ( r/10 4 K)~ 1/2 is the C ase B 
recombination rate of hydrogen (Maselli et al. 20031. Nor- 



Mrec = (47r/3)a' 2 'n| f0 r 3 ir , and adopting a dimensionless 



time variable r 



t / ire 



a 2 'riHot, equation (|9l) can be 



written in a dimensionless form, 



dx 
d^ 



Q(t) 



2>Nre C r)(x)x 2 



-xr](x). 



(10) 



Equation ( 10 I describes the expansion of the H n region into 



the stratified gas distribution within the halo. As the density 
decreases outwards, the IF accelerates and eventually exits 
the virial radius at a time t out = t(x — 1). Until t ou t all 
the ionizing photons are absorbed inside the halcP] hence, 
f e3C « 0. However, for t > t out , ionizing photons will be only 
used to keep the halo ionized by balancing recombinations 
within it, but a large fraction of them will be able to escape 
into the IGM, thus making f esc > 0. We then express the 
net escape fraction as follows: 



fesc — 



J °°drQ(t) ■ 



(11) 



As we will see shortly, though, either the IF is efficiently 
confined by the density, dx/dr — > 0, or the blow out will 
occur on a very short time scale thus making f eS c w 1. We 
will refer to these two different situations as a "confined" or 
"unconfined" IF in the rest of the paper. 



3 RESULTS 

We solved numerically equation ( |10[ ) for a number of halo 
masses with V — 1 in the range M = 10 s -10 11 Mq at two 
selected redshifts, z = 9 and z = 14. Figures [I] and [2] show 
the results for /* = 0.2%, a natural value for for nearly 
primordial star- forming halos ( Barkana fc Loeb|2001|[Ciardi 
fc Ferrara|2005 l. 

IFs in these early halos evolve rapidly to the strong ra- 
diative input by the stars and reach a radius ri = r V i r within 
~ 1 Myr, yielding / esc = 0.98 — 0.90 for halo masses up to 
Mr;2x 10 10 M , beyond which the IF is confined within 
the halo, and hence / esc = 0. It requires fine tuning for a 
single object to achieve intermediate / e3C values; rather, the 
process favors on-off states implying that the escape fraction 
of photons from galaxies can be intermittent: depending on 
variations of /* in different bursts (say, due to the increasing 
gas metallicity or decreasing gas fraction due to mass loss 
from supernova-driven winds) the same halo might switch 
from being dark to bright above the Lyman limit. Note also 
that for the largest IF-confined halos, the H n region starts 
to shrink at t Si 3 Myr as a result of recombinations and the 
decline in the supply rate of ionizing photons. 

The situation is dramatically different at z = 14 (for 
the same /*): at that epoch all halos confine their IF and 
have fesc = 0. This is due to the increase in the gas den- 
sity and the larger concentration factor at a fixed mass. A 
more global view of the dependence of f esc on halo mass 
and star formation efficiency at z = 9 is given by Fig. [3] 
The abrupt transition from f esc ~ 1 to zero is very weakly 
dependent on halo mass, except for the narrow range of 
—2.8 < log/* < —2.2 in which larger halos tend to be more 



malizing by an effective recombination rate in the halo, 



1 http:/ /www. stsci.edu/science/starburst99/ 



2 Here we ignore the possibility of a highly inhomogeneous 
medium in which low density channels for escape exist. 
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Figure 1. Time evolution of the ionization front radius, rj, at 
z = 9 within halos of different total mass 

with j = 0, .., 9 from the top curve (M = 10 s M Q ) to the bottom 
one (M = 10 Mq). The assumed star formation efficiency is 
/* = 0.002 and the gas temperature factor V = (T^r/T) 1 / 2 = 1. 
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Figure 2. Same as Fig.[T] but for z = 14. 



opaque to ionizing photons. Above this range, f esc — 1 in- 
dependently of halo mass and star formation efficiency, and 
below this range, f esc — 0. 

Figure [4] shows the examples with V > 1 introduced 
in Table 1, for which T = 2 x 10 4 K, corresponding to 
V 2 = (1.04,4.82,22.39, 103.96) for halos of mass M/M = 
(10 8 , 10 9 , 10 10 , 10 11 ) and hence resulting in considerably 
denser and more compact gas configurations. The IF front 
is now largely confined within all halos with M > 10 s Mq 
even if the star formation efficiency has been increased by a 
factor of 50 with respect to the case with V = 1 shown in 
Fig.[T] In addition, the IF escapes the virial radius in a very 
short tim^] t out ~ 10 4 yr. 

The different behavior relative to the case of V = 1 is 



We have checked that the IF speed is subluminal, drj/dt < c 
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Figure 3. Map of f esc values as a function of halo mass and star 
formation efficiency at z = 9 assuming V = (T v i r /T) 1 / 2 = 1. The 
color scheme defined by the vertical bar on the right hand side. 



caused by the associated high gas density (see Table 1). The 
small scale height of the gas distribution, correspondingly 
reduced by a analogous factor oc V 2 , forces the steep ac- 
celeration seen just prior to blow out. In analogy with the 
V = 1 case, we show the dependence of / esc on halo mass 
and star formation efficiency at z = 9 (Fig. ^ . In this case 
the / e3C = 1 area is much narrower and limited to small 
halos with M < 10 8 ' 7 M Q and /* > 10~ 3 . Thus, according 
to our model, it is only the smallest halos that are truly 
able to inject ionizing photons in the IGM and therefore 
substantially contribute to reionization. Larger halos are in- 
stead able to confine their H n regions very effectively. At 
z — 14 the / esc = 1 region is qualitatively the same, but its 
area shrinks because of the larger background density and 
the increased concentration parameter at a fixed halo mass. 

The existence of a minimum value of /* w 10~ 3 to allow 
the escape of photons can be understood based on a simple 
argument. This minimum number of ionizing photons pro- 
vided to each hydrogen atom in the halo must balance the 
number of recombinations it undergoes during the lifetime 
of the starburst, approximately equal to to — 10 6 ' 6 yr (see 
equation |8J. This can be expressed by the following inequal- 
ity: 



/* > TT—Z—( 1S ' K )^bPc(z) 



N~f /imp 



0.7 x 1Q" A — 



5 x 10 3 



10 



(12) 



(13) 



where C = (n H ) / (tih) ~ 3 is the clumping factor account- 
ing for the gas density structure inside a NFW halo. 

The results for halos with M < 10 8 M in which 
T mr < 2 x TO 4 K at z = 9 need some extra attention. In 
fact, by fixing the gas temperature to 2 x 10 4 K, the gas 
might have never been accreted in the first place; this is the 
case if for example the halo is located within an already ion- 
ized region. To correct this problem, we have set V = 1 in 
these mini-halos; of course, this requires that some cooling 
agent, such as H2 , HD or LiH, would cool the gas in these 
systems to make stars. 
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Table 2. Values of various quantities related to the global evolution of the mass-averaged escape fraction. See text for more details. 
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4.852e+08 
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1.169c+07 


1.169e+07 
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1.130e-01 


5.936e-02 


4.747e-01 


3.000e-02 


z 


= 9 


1.169e+07 


1.169e+07 


5.818e+08 


1.130e-01 


4.374e-02 


6.130e-01 


1.000e+00 


z 


= 9 


1.169e+07 
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1.130e-01 


1.130e-01 


0.000e+00 


1.000e-04 


z 


= 14 
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5.373e-02 


1.471e-02 


7.263e-01 


3.000e-02 
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Figure 4. Same as Fig. ^ but for halos with mass in the range 
M = 10 8 +2j '/ 30 M©, with j = 0, ..,9 from the top curve (M = 
10 s M ) to the bottom one (M = 3x 10 s M Q ). The star formation 
efficiency is /* = 0.1 and the gas temperature has been fixed to 
T = 2 X 10 4 K, giving a variable V temperature factor (see text 
for details). 



3.1 Global redshift evolution 

Next we would like to constrain the redshift evolution of 
f eac (z) when averaged over the entire galaxy population. 
The dark matter halo mass function, n(M, z), is well 



described by the Sheth & Tormen ( 2002 1 form 



n(M, z)dM = A 1 + 



..llq 



, - r -ttt exp dM, 

n M dM v \ 2 J 

(14) 

with A = 0.322, q = 0.3, v' = sfav, a = 0.707, v = 
S c / D(z)a(M), with S c = 1.686, D(z) being the linear growth 
factor and a(M) being the r.m.s. amplitude of density fluc- 
tuations on a mass scale M. 

We define the mass-averaged escape fraction, (fesc) at 




7.0 7,5 8.0 8,5 
loRfAf/Af R l 

Figure 5. Same as Fig. [3] but for a fixed gas temperature T 
2 X 10 4 K. 



(15) 



a given /* and redshift z as follows: 

(f ){z f ) _ x F[>M f0 {z,U)] 

where F(> PI, z) = /~ nMdM/ nMdM is the collapse 
fraction of all matter in haloes with mass exceeding M at 
redshift z. 

Equation ( |15[ ) contains two important characteristic 
masses. The first, Mfo(z, /*), is defined as the halo mass 
above which, for a given redshift z and star formation ef- 
ficiency /*, no ionizing photons can escape, i.e. f esc = 0. 
The results presented above allow us to precisely determine 

M/o (*,/*)• 

The second mass scale, M s f{z), denotes the redshift- 
dependent minimum mass of star-forming haloes which is set 
by radiative feedback. The evolution of M s f is determined 
by two distinct radiative feedback processes. The first one 
is related to the increase of the cosmological Jeans mass in 
progressively ionized cosmic regions; as a result, the infall of 
gas in haloes below a given virial temperature is quenched. 
The evolution of M s f(z) depends on the details of the reion- 
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ization history (Gncdin 2000; Okamoto et al. 2008; Schnci- 
der et al.|2008 1. Guided by these findings, we adopt the the 
value M 3f = M(T vlr = T 2 ) with T 2 = 3 x 10 4 K after the 
end of reionization, here assumed to occur at z re i — 6. 

Before reionization (at least in the mostly neutral re- 
gions) a second type of feedback, involving the photodisso- 
ciation of H2 molecules by the Lyman- Werner (LW) back- 
ground photons, becomes important. As H2 is the primary 
coolant for minihalos (T V i r < 10 4 K), fragmentation of the 
gas to stars could be suppressed in these objects through 
H2 dissociation by the UV background ( |Haiman" et al.|1996| 



Ciardi et al. 2000; Kitayama et al. 



2000| IM achacck et al.| 
" 2004| ), 



2001[ ). According to |Dijkstra et al.| (|2004), z ~ 10 mini- 
halos with T vir = Ti » 2x 10 J K can self-shield H2 and cool, 
hence we use this value as the lower threshold for star forma- 
tion. During reionization the interplay between the photo- 
ionization heating and LW-feedback is complicated (see the 



discussion in Mesinger & Dijkstra ( 2008 1 and Fig. 25 of Cia- 



rdi & Ferrara ( 2005 1 ) and and the evolution of M 3 f during 
this epochs is uncertain. To circumvent this uncertainty, we 
follow the phenomenological approac h of |Salvadori fc Fer-| 



(20091 and ISalvadori & Ferrara (20121 who suggested 



that the metallicity distribution and iron-luminosity rela- 
tion of the Ultra Faint dSphs in the Milky Way halo can be 
explained if these objects are relic mini-halos that formed 
their stars prior to reionization. With this perspective, one 
can reconstruct M s f(z) by fitting local data. This approach 
yields the functional form, 



M sJ {z) = M[T viT = (T 2 - 7\)e ( 



+ Ti] 



(16) 



with Az = 4, which we adopt here. 

In general M/o(z, /*) > M s f(z). However, for suffi- 
ciently low values of /*, Mfo(z, /*) = M s f(z) and (/ esc ) = 0. 
To illustrate this point (see Table we examine the results 
at z — 9 in Figure [5] For /* = 0.03, stars are formed only 
in halos more massive than M s f = 1.17 x 10 7 Mq. Ioniz- 
ing photons can only escape from halos less massive than 
Mfo = 2.19 x 10 8 Mq. As a result, there is a narrow range 
of objects that contribute to the ionizing photons budget 
in the IGM. Below that range stars cannot form and above 
it the gas is sufficiently dense to confine the H 11 region. 
When weighted with the halo mass function, though, the 
collapsed mass fraction above (M s f , Mfo) is (0.113,0.059), 
yielding the large value (f esc ) (z = 9, /* = 0.01) = 0.847. 
Since H 11 regions in abundant low mass halos are unbound, 
the average escape fraction is high. The width of the range 
progressively shrinks as star formation becomes less efficient. 
This is evident from Table [2] where for 2 = 9 and /* = 10~ 4 , 
M/o(z,/*) = M s f(z) and the escape fraction drops to zero. 
These results underline the importance of the possible evolu- 
tion in the star formation efficiency, which we address below. 




Figure 6. Redshift evolution of the mass-averaged escape frac- 
tion for different assumptions concerning the star formation effi- 
ciency. Solid lines refer to the redshift-independent values /* = 
0.03, 0.003 as indicated by labels; the shaded area refers to /*(z) 
(shown by the shaded cyan area and multiplied by a factor 100) 
derived from the stellar mass density evolution by Gonzalez et al. 
| |2011[ l based on equation \17\ . 
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Figure 7. Cumulative fraction of ionizing photons produced by 
halos below a given halo mass M. Lines refer to a constant 
star formation efficiency /* = 0.003; shaded regions refer to 
the upper limit of the redshift-dependent /* (equation |17| based 
on |Gonzalez et al.] ( |2011[ l SMDs. For both sets of curves, solid 
(dashed) lines refer to 2 = 8 (z = 10); for /* = 0.003 the case for 
z = 12 (dotted) is also shown. 



The redshift evolution of the results can be understood 
in terms of a balance between two opposite trends. On the 
one hand, the halo central densities tend to become larger 
with increasing redshift, providing a more efficient confine- 
ment of the IF and acting to decrease (/ eS c). On the other 
hand, more dwarf galaxies and mini-halos are able to pro- 
duce ionizing photons as the critical mass M s f decreases 
with increasing z. The overall evolution for the two spe- 
cific values /* = (0.03, 0.003) can be seen in Figure [6] (solid 
curves). For the high case, the escape fraction continues 



to increase with redshift, reaching (fesc) = 0.83 at z — 20. 
If instead /* — 0.003, the escape fraction reaches a peak 
value of 0.59 at z = 11.5 where it starts to decrease until it 
reaches zero at z = 17.5 when the density effect dominates 
over the enhanced photon production rate. Interestingly, a 
nearly redshift-independent /* above a threshold of w 0.1% 
(see Fig. [5| always leads to d(f esc )/dz > during the final 
phase of reionization. Although a rather weak mass depen- 
dence of /* for galaxies of stellar mass M* > 3 x 10 7 Mq is 
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inferred from local SDSS data (Dekel & Woo 20031, there 
are indications of a decreasing star formation efficiency for 
high-z galaxies. Recently, different groups ( Stark et al.|2009 



Labbe et al.|2010 Gonzalez et al.|[201 1 1 have obtained stel 



lar masses of galaxies in the redshift range 4 < z < 7 
from spectral energy distribution fitting to rest-frame op- 
tical and UV fluxes from Hubble-WFC3/IR camera obser- 
vations of the Early Release Science field combined with the 
deep GOODS-S Spitzer/IRAC data. From this data it has 
been possible to reconstruct (see Fig. 4 of Gonzalez et al. 



( 2011 1) the evolution of the comoving Stellar Mass Density 
(SMD), which we define here as p*(z). The main result is 
that p*(z = 6) = 10 7 M Q Mpc" 3 ; in addition, the SMD 
grows with decreasing redshift as oc 

(1 + z) -3.4±0.8 We 

can 

exploit this observational result to obtain a robust estimate 
for f*(z) based on the relation, 



/*(*) 



f b F[>M sf (z)] p c 



(17) 



In other words, we tune the star formation efficiency so that 
baryons in collapsed halos account for the abundance of stars 
observed at z < 7. Note that the |Gonzalez et al.| ( [20lT| SMD 
has also been corrected for incompleteness at low masses. 
The range of allowed values of /* implied by the data and 



obtained through equation ( 17 1 are plotted as a shaded area 
in Fig. [6] The data suggest /* = 0.003 at z = 6 (for display 
reasons we have multiplied the values by a factor 100 in 
the Figure) decreasing thereafter by a factor of about three 
above z = 10. 

The main consequence of a suppressed star formation 
efficiency in low mass objects is to introduce a new interme- 
diate mass scale M s f < M/i < M/o such that f esc > only 
in the interval [M/i,M/o]- Stated differently, this is an evi- 
dence of the fact that even objects that are above M s f and 
can form stars, are now more efficient at confining their IFs; 
therefore their escape fraction becomes zero. Table [2] shows 
how Mfi gradually approaches M/o towards high redshifts, 
closing the non- vanishing escape fraction window. 

Another interesting question is how the global produc- 
tion of ionizing photons is distributed among halos of differ- 
ent masses. We quantify this through the cumulative fraction 
of ionizing photons produced by halos below a given mass, 

r-M 

/ 7 (<M,z)cx/ f eac (M',zJ,)f ir (M af ,z)M'dM', (18) 



as shown in Figure [7] Due to their large escape fraction, ha- 
los with masses M 10 8 Mq are the dominant contributors 
to the cosmic ionizing photon budget. This holds for the 
more realistic redshift-dependent /*, but also for the con- 
stant /* case. Thus, the low-mass galaxies and mini-halos 
might have dominated reionization. In particular, at z — 10 
the mass corresponding to T = 10 4 K limiting the mini-halo 
mass range from above, is ~ 3 x 10 7 Mq. Figure [7] shows that 
independently of the star formation efficiency prescription, 
we predict that mini-halos contribute ~ 40% of the total 
emissivity. 

Reionization models suggest an increasing emissivity to- 
wards hi gh redshifts (|Bolton fc Haehnelt||2007| |Mitra et aT] 
2012a|b|[Kuhlen fc Faucher-Gi guere|20l"2"| , whose most nat- 
ural explanation is in terms of a / esc increasing with redshift 
as we find here. For example, |Mitra et al.| ([2012b) find that: 
(i) the escape fraction increases from (f esc ) = 0.068tg;S!g 



at z=6 to 0.179±q ?|a 

(fesc) > 0.146. 



at z = 8; and (ii) at z 



10 



4 DISCUSSION 

We have calculated the escape fraction of ionizing photons 
from starburst galaxies during reionization, assuming that 
galactic gas cools down to a temperature floor (~ 10 4 K for 
atomic hydrogen and ~ 2 x lO^K for self-shielded molec- 
ular hydrogen). We have found that most of the escaping 
ionizing photons originates in the lowest mass galaxies with 
Tvir ~ 10 4 K. Due to the shallowness of the gravitational- 
potential well in these galactic halos, the gas maintains 
a low density and its recombination rate is low. This al- 
lows the ionization front to break through the virial radius 
rapidly, resulting in a high escape fraction of the emitted ion- 
izing radiation. As reionization progresses, the intergalactic 
medium is photo-heated and the abundance of galaxies with 
T v i r ~ 10 4 K declines, leading to a decrease in (/ eS c). This 
physical mechanism may account for the increase of (/esc) 
with redshift inferred by Mitra et al. (2012b I . 



Our study is consistent with previous studies of low 



mass galaxies by Wise & Cen (20091, who found a high 



escape fraction (> 0.1). Low mass galaxies do not form a 



disk (Pawlik et al. 20111, which could trap effectively ion- 
izing photons. The relative importance of minihalos for the 
global production rate of ionizing photons is determined by 
the product f eac x /*. Here we have used a phenomenologi- 
cal calibration of /* . In order to produce the observed stellar 
density at z = 6-7 by minihalos, /* has to be £ 0.001. This 
value lies in the middle range resulting from |Wise fc Cen| 
( 2009 ) idealized halo simulations, but they find the star for- 
mation efficiency in halos with M < 10 7 Mq to be smaller. 
We find that the required increase of / esc with redshift can- 
not be achieved if minihalos are sterile or form stars ineffi- 
ciently. 

Our simplified calculation ignored angular momentum 
which drives cold gas to a rotating disk configuration. Our 
model allowed cold gas to condense to arbitrarily small scales 
without making a disk, and so one might worry that it under- 
estimates / eBC . However, previous calculations have demon- 
strated that on the scale of galactic disks the escape fraction 



is very small (Wood & Loeb||2000 


Gnedin et al. [12008 Ra- 


zoumov & Sommer-Larsen 2006 


Wise & Cen 20091, and 



so our inclusion of angular momentum could not increase 
/esc significantly. In particular, the conclusion that / eS c is 
high in low-mass halos must be robust. Although our sim- 
plified model predicts zero escape fractions from brief star- 
bursts in massive halos, the actual value of / eS c would be 
finite in reality due to persistent star formation histories 
and winds driven by supernova feedback or Lya radiation 
pressure ( Dijkstra fc Loeb|2008 l. Our results are in qualita- 



tive agreement with Gnedin et al. ( 2008 \ and Razoumov & 



Sommer-Larsen ( 2006 \ who simulated more massive galax- 
ies (M ^ 10 iu Mq) and inferred an extremely small escape 
fraction values f esc ;$ 1%), even in situations where the star 
formation rate is high (~ 10 Mq yr~ ). This results from 
the high gas densities in massive galaxies due the presence 
of a rotationally-supported, geometrically-thin disk. 

Finally, we note that the trend we find for the increase 
of the escape fraction with redshift should be enhanced if 
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the initial mass function of stars was tilted towards more 
massive stars. This is because a higher production rate of 
ionizing photons would allow the ionization front to ionize 
the surrounding hydrogen blanket even faster than we cal- 
culated. 
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